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Abstract 

By applying the multifractal detrended fluctuation analysis to the high-frequency 
tick-by-tick data from Deutsche Borse both in the price and in the time domains, we 
investigate multifractal properties of the time series of logarithmic price increments 
and inter-trade intervals of time. We show that both quantities reveal multiscaling 
and that this result holds across different stocks. The origin of the multifractal char- 
acter of the corresponding dynamics is, among others, the long-range correlations in 
price increments and in inter-trade time intervals as well as the non-Gaussian distri- 
butions of the fluctuations. Since the transaction-to-transaction price increments do 
not strongly depend on or are almost independent of the inter-trade waiting times, 
both can be sources of the observed multifractal behaviour of the fixed-delay returns 
and volatility. The results presented also allow one to evaluate the applicability of 
the Multifractal Model of Asset Returns in the case of tick-by-tick data. 
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1 Introduction 



From the perspective of the classical financial market models, behaviour of 
the consecutive price and index fiuctuations does not present any significant 
time autocorrelations except for short time scales up to several minutes. As 
regards the dynamical character of these fiuctuations, they are considered 
as being to a good approximation the fractional Gaussian noise [1], with very 
small and negligible probability of the occurrence of non-Gaussian large jumps 
in the index or the share price. As a natural consequence, the stock mar- 
ket data is expected to present only monofractal properties. However, these 
widely-used models do not describe the processes underlying the evolution of 
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financial data with satisfactory precision. The so-cahed financial stylized facts 
comprising, among others, the non-negligibic fat tails of log-return distribu- 
tions, volatility clustering and its long-time correlations, anomalous diffusion 
etc. [2,3,4,5] counter the above-mentioned fundamental assumptions of market 
dynamics challenging their applicability in practice. That the financial dynam- 
ics is more complex than it is commonly assumed can also be inferred from 
a number of recently-published papers discovering and exploring the multi- 
fractal characteristics of data from the stock markets [6,7,8,9,10], the forex 
markets [7,10,11,12,13] and the commodity ones [14]. The concept of multi- 
fractality was developed in order to describe the scaling properties of singular 
measures and functions which exhibit the presence of various distinct scal- 
ing exponents in their different parts [19,20]. Soon the related formalism was 
successfully applied to characterize empirical data in many distant fields like 
turbulence [21,22], earth science [23], genetics [24,25,26], physiology [27,28,29] 
and, as already mentioned, in finance. The problem of detecting multifractality 
in real data is delicate, however. There are models based on fractal processes 
which are able to mimic the real multifractal evolution of markets being either 
multifractal or monofractal [15,18,30,31,32]. Moreover, as it has been pointed 
out ([18,33]), the power of commonly-used tests of multifractality is limited, 
because they cannot effectively distinguish between the two types of fractal 
behaviour of the financial (but perhaps also other) data. One important source 
of this difficulty is the presence of non-Gaussian tails in the distributions of 
data (e.g. truncated Levy [33]), the fact which is ubiquitous in finance. Thus, 
all conclusions drawn from multifractal analysis have to be interpreted with 
care. 

In the present paper we analyze data from the German stock market focusing 
on their fractal properties. We apply the multifractal detrended fiuctuation 
analysis which is a well-established method of detecting scaling behaviour 
of signals. By exploiting the character of the high-frequency transaction-by- 
transaction recordings for the most liquid stocks belonging to DAX, we are 
able to inspect not only the properties of share price fluctuations, but also the 
properties of time intervals between consecutive trades (waiting times). The 
majority of analyses carried out so far was devoted to time series of the returns 
calculated with some fixed time delay At. According to the Multifractal Model 
of Asset Returns introduced by Mandelbrot and others [15,16,17,18,31], the 
source of multifractality in the returns is a deformation of time 0{t), which 
takes place due to the fact that at the microscale the so-called business time 
is "paced" by transactions rather than any constant time units. In this model 
the increments of so-deformed time can be directly related [4,31] to the fixed- 
At volatility which depends both on the trading activity (i.e. ~ ^ N^^^\t) , 
where A^(^*) stands for the number of transactions in At) and on the variance 
W^'^^^t) of the price change in the individual transactions over At. Looking 
deeper into the microstructure of the market, as these quantities are equal to. 
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respectively, the reciprocal of the average waiting time in At and to the average 
squared price increment in At, one may ask whether, at the microscopic level, 
the multifractal properties of 6{t) originate from the price behaviour, from 
the waiting times fluctuations or from both of them. This is our motivation 
for performing the present analysis. We find the German stocks particularly 
suitable for such a study due to the fact that the moments of transactions 
are recorded with high precision (0.01 s), which almost completely eliminates 
falsely simultaneous transactions. Such an analysis for the liquid stocks from 
NYSE, although desired, cannot be so successfully carried out as the time 
resolution of the recordings is poor (Is) and too many transactions are forced 
to be simultaneous, distorting the underlying real dynamics. 



2 Methods and data 

There are two possible procedures of analyzing multifractal properties of a 
time series. The first one uses the continuous wavelet transform and extracts 
scaling exponents from the wavelet transform amplitudes over all scales [34]. 
This method is more computationally demanding and, as our preliminary tests 
showed, the stabihty of results for our data is not satisfying. Therefore, for the 
present study we prefer to employ the multifractal version of the detrended 
fiuctuation analysis method (MF-DFA) [24,35]. 

Given the time series of price values Ps{ts{i)),i = 1, A^s of a stock s recorded 
at the discrete transaction moments ts{i), one may consider two independent 
random processes defined by price Ps{i) and time ts{i) or, alternatively, loga- 
rithmic price increments ^ gsi'i) = lii(ps(i + 1)) —ln(pg(i)) and time increments 
(waiting times) Ats{i) — ts{i + 1) — ts{i)- 

For the time series of the log-price increments Gg := {gs{i)}i=i,..,Ns (and analo- 
gously for the waiting-times series Tg :— {ts{i)}i=i,...,Ns): oiie needs to estimate 
the signal profile 

Y{i) = E {gs{k)- <gs>),i^l, (1) 
fe=i 

where < ... > denotes the mean of Gs- V(i) is divided into Ms disjoint 
segments of length n starting from the beginning of Gg- For each segment 



Throughout this paper wc intend to use the expression "price increments" in- 
stead of "returns" to underline the fact that we study the price differences between 
consecutive trades which occur irregularly in time, while in literature the "returns" 
are usually associated with a constant time step. 
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— 1, Ms, the local trend is to be calculated by least-squares fitting the 
polynomial Pj'^ of order / to the data, and then the variance 

1 

F\v, n) = - - l)n + j] - pI%j)Y . (2) 

In order to avoid neglecting data points at the end of Gs which do not fall into 
any of the segments, the same as above is repeated for Mg segments starting 
from the end of Gg (i.e. finally one has 2Afs segments total and the same 
number of F^'s). The polynomial order / can be equal to 1 (DFAl), 2 (DFA2), 
etc. The variances (2) have to be averaged over all the segments v and finally 
one gets the gth order fluctuation function 

( 1 2Ms . l/q 

In order to determine the dependence of Fq on n, the function Fq{n) has to 
be calculated for many different segments of lengths n. 

If the analyzed signal develops fractal properties, the fluctuation function re- 
veals power-law scaling 

Fq{n) - n^^'i^ (4) 

for large n. The family of the scaling exponents h{q) can be then obtained by 
observing the slope of log-log plots of Fq vs. n. h{q) can be considered as a gen- 
erahzation of the Hurst exponent H with the equivalence H = h{2). Now the 
distinction between monofractal and multifractal signals can be performed: if 
h{q) = H for all q, then the signal under study is monofractal; it is multifrac- 
tal otherwise. By the procedure, h{q),q < describe the scaling properties 
of small fluctuations in the time series, while the large ones correspond to 
h{Q.)^ ? > 0. It also holds that h{q) is a decreasing function of q. 

By knowing the spectrum of the generalized Hurst exponents, one can calcu- 
late the singularity strength a and the singularity spectrum /(a) using the 
following relations (e.g. [35]): 

a^h{q)+qh'{q) f{a)^q[a-h{q)] + l, (5) 

where h'{q) stands for the derivative of h{q) with respect to q. 

Our analysis was performed on the time series of tick- by-tick recordings for 
the 30 DAX stocks comprising the over-two-years-long interval between Nov 
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Fig. 1. Time course of time series of logarithmic price increments gs{i) (upper panel) 
and time series of inter-trade waiting times Ats{i) (lower panel) for an exemplary 
stock (BMW). 

28, 1997 and Dec 31, 1999. The time series were approx. 250,000 points long 
on average (almost 500 transactions daily) with the minimal length of 63,000 
(Karstadt Quelle, KAR) and the maximal one of 588,000 (Daimler-Chrysler, 
DCX). Firstly, we removed all the overnight price increments, because they 
correspond to extremely long inter-trade intervals and largcr-than-usual price 
jumps (such effects were identified to introduce distortion of the financial 
scaling [36]); we do not remove zero increments from Gg and zero waiting 
times from because the corresponding zero-intervals arc so short that tlicy 
do not alter the results. The so-preprocessed data was subject to the MF-DFA 
analysis with the polynomials P^^^ (MF-DFA(3)), owing to the fact that for 
the present data they optimally extract the fluctuations whose scaling is to be 
analyzed (Eq. 4). 



3 Results 



Figure 1 shows the time course of Gg (upper panel) and Tg (lower panel) for the 
complete interval under study. There is a clear difference between the series. 
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Fig. 2. Log-log plots of the q-ih. order fluctuation Fq for time series of price incre- 
ments as a function of segment size n for different values of q between -10 (bottom 
line) and 10 (top line). Six stocks are shown ordered according to the slope spread 
between the extreme values of q (left to right and top to bottom). In each panel 
only a region of the scaling regime is shown and a non-scaling part for small n has 
been cut off. The largest value of n depends on the time series length being different 
for each stock. Scaling regions allow one to estimate h{q) according to Eq. 4. 

because while the distribution of the former is symmetric around zero, the 
latter cannot assume negative values and, hence, its distribution is skewed. 

First we shall present results of our study of the price increments data. Fol- 
lowing Eq. (4), for each stock we created double- logarithmic plots of Fq{n) 
for the various segment lengths in the range 30 < n < Ns/h data points and 
for various choices of q. It is widely assumed that the amplitude of the price 
fluctuations scale according to the inverse cubic power law, which in principle 
implies infinite moments for g > 3. In real data the large-order moments are 
therefore inevitably affected by the finite-size effects. However, we decided to 
include the values of q as high as 10 purely to improve readability of the plots 
of the multifractal and singularity spectra and to make any comparison of the 
spectra's widths easier; this step does not reduce the significance of our results. 
The plots of Fq{n) for six representative stocks from different market sectors 
are collected in Fig. 2 and ordered according to the decreasing slope spread 
between g = — 10 and 10 (with the step of 0.2 for small |g|'s and of 0.5 for 
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Fig. 3. Multifractal spectra expressed by r(g) = qh{q) — 1 for price increments; 
dot-dashed line in each panel presents a monofractal Gaussian noise regime. A 
nonlinear behaviour of r(g) can be considered a manifestation of multiscaling. The 
same stocks as in Fig. 2 are shown. 



larger ones). In order to indicate the range of n's used for fitting the exponents 
we present the scaling regions only and cut off the regions characterized 
by the lack of scahng (this happens for small n's, probably due to the existence 
of long periods of constant price of a stock). The largest difference between the 
slopes is observed for Adidas Solomon (ADS, top left panel) and the smallest 
one for BMW (bottom right panel). The plots for different companies show 
noticeable differences in the range of n with the scaling behaviour; the widest 
range of over three decades is observed for moderate values of q for ALV (Al- 
lianz) and SIE (Siemens). In each case, for extremely small and extremely large 
n's the scaling breaks down as expected from statistical considerations [35]. 
In order to better visualize the scaling character of the data, in Fig. 3, we 
show the corresponding multifractal spectra. Instead of h{q), we display its 
function r(g), defined by the relation r(g) = qh{q) — 1. Monofractal signals 
{h{q) =const) are associated with a linear plot T(g), while multifractal ones 
possess the spectra nonlinear in q. Keeping in mind the already-mentioned 
limitations of the method ([17]), our calculations indicate that the time se- 
ries of price increments for all companies can be of the multifractal nature. 
Consistently with the log-log plots on the left, the highest nonhnearity of the 
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Fig. 4. Singularity spectra / (a) vs. a for price increments and for the same stocks as 
in Fig. 2. The maximum value of each spectrum gives the most common singularity 
strength ao- The vertical dotted lines denote q = —3.0 (right) and q = 3.0 (left). 

spectrum and the strongest multifractality are attributes of ADS and DBK 
(Deutsche Bank), and the smallest nonlinearity and the weakest multifractal 
character correspond to BMW and BAY (Bayer). The nonlinearity of r{q) is 
confined to the central range of g's around q — and for larger values of |g| 
the behaviour of T{q) is almost linear (due to the finite size of the sample). 

The multifractal nature of the data can also be expressed in a different manner, 
i.e. by plotting the singularity spectra f{a) (Eq. 5). It is a more plausible 
method because here one can easily assess the variety of scaling behaviour 
in the data. Figure 4 displays such spectra for the same stocks as in Figs. 2 
and 3, with their presentation order being preserved. As before, the richest 
multifractality (the widest f{a) curve) is visible for ADS (Aa := ag=_io — 
aq=io ~ 0.5) and DBK (Act ~ 0.45), the poorest one for BAY and BMW 
(Act ~ 0.3). The maxima of f{a) are typically placed in a close vicinity of 
a = 0.5 indicating no significant autocorrelations exist. It is worth mentioning, 
that the singularity spectra for the price changes between two consecutive 
transactions closely resemble their counterparts for the fixed-At log-returns. 
We do not show the corresponding results graphically, but the widths of the 
/(a) spectra are similar in both cases. 
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Fig. 5. Comparison of the original and randomized price increments: original (solid), 
reshuffled (dashed) and surrogate (dotted) time series. Multifractal spectra r(g) (left 
column) and singularity spectra f{a) (right column) are presented for three different 
stocks. 

The multifractal character of price fluctuations can originate from the exis- 
tence of the long-range correlations in the price increments (via volatility) as 
well as from their non-Gaussian distributions [33]. The possible influence of 
each of these factors can be detected by a proper modiflcation of the data. 
The long-range autocorrelations can be completely erased by randomly reshuf- 
fling the original time series and the non-Gaussianity of the distributions can 
be weakened by creating the phase-randomized surrogates [37]. In the latter 
case we exploit the fact that the price increments distributions are unstable 
in the sense of Levy, which leads to their convergence to a Gaussian under 
the discrete Fourier transforms. Figure 5 shows three examples of r{q) (left 
column) and f{a) (right column) spectra for the original (solid), reshuffled 
(dot-dashed) and phase-randomized (dashed) data. Both the widths of the 
f{a) spectra in each case are much smaller and the nonlinearity of T(g)'s is 
much weaker for the modified signals than for the original ones. This behaviour 
of the reshuffled signals conflrms that the persistent autocorrelations play an 
important role in multiscaling of the price increments. However, the spectra 
for the surrogates are typically much narrower than for the reshuffled data 
which can be interpreted as an evidence of the influence of extremely large 
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Fig. 6. Log-log plots of Fq for waiting times as a function of n (— 10 < g < 10). This 
figure is the analog of Fig. 2 for Gs] the same six stocks are shown. 

non-Gaussian events on the fractal properties of the signals. Keeping in mind 
the difficulties in precisely calculating the scaling exponents h{q) for finite 
time series, one may interpret the narrow curves for the modified signals in 
Fig. 5 as the manifestation of their relatively, but not precisely, monofractal 
character. 

To this end, we concentrated on multifractal properties of the time series of 
logarithmic price increments. Our results go in parallel with earlier analyses 
of other groups, which managed to show multifractality in the stock market 
data [8,12,14]. Once again, we stress that what distinguishes our approach is 
that we did not sample the data with fixed time interval, but rather look into 
the tick-by-tick data and constructed the time series of price increments for 
variable inter-trade time intervals (thus unfolding time). Now we shall go back 
to the original time axis, and study properties of the inter-trade time intervals 
forming the series Tg. 

Figures 6, 7, and 8 are fully analogous to Figs. 2, 3, and 4, respectively, but 
for the waiting times instead of the price increments. As before, we observe 
significant difference between the slopes for the most negative and the most 
positive g's with this strength being different among the stocks. The scahng in 
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Fig. 7. Multifractal spectra r(g) = qh{q) — l for waiting times. A nonlinear behaviour 
of T{q) can be considered a manifestation of multiscaling. 

Fig. 6 is good for small and moderate values of |g|, but is poorer for large |g|'s 
for the majority of companies (the effect of finite size and noise). The r(g) 
spectra in Fig. 7 show strong nonlinearity, and although it is confined to middle 
range of g's, this nonlinearity is stronger than for the price increments. By 
comparing Figs. 2-4 and 6-8 we see that there is no systematic evidence of the 
relations between the properties of Gs and for the given companies. Fig. 8 is 
a counterpart of Fig. 4 and presents the singularity spectra for the time series of 
the waiting times. All the spectra can be interpreted as multifractal. The main 
differences are the larger widths of /(q;)'s, their asymmetry and the positions 
of maxima at approx. a > 0.8 instead of 0.5 for price fluctuations. The widest 
spectra correspond to ADS and ALV, the significantly narrower ones to BAY, 
in different order than for Gg- The shift of the maxima can be related to 
long-range correlations between consecutive waiting times, leading to a strong 
persistence (see also Fig. 9). However, despite the fact that both the spectra for 
Gg and for Ts are multifractal, they cannot be directly compared to each other. 
The crucial factor here is that they represent processes of different character: 
signed and unsigned, respectively. There exists a method of rcscaling the /(«) 
spectra of the unsigned (or signed) process (as it was proved in [16]), but 
it requires that the underlying processes X{t),Y{t) are strongly interrelated 
through the fractional Brownian motion X{t) — BH[Y{t)] with some Hurst 
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Fig. 8. Singularity spectra /(a) for waiting times; the maximum value of each spec- 
trum gives the most common singularity strength ao- The vertical dotted lines 
denote q = -3.0 (right) and q = 3.0 (left). 

exponent H (then simply fx{c() = fyia/H)). However, in the context of our 
data, this assumption seems to be violated: Gs and Ts can in principle be 
independent or, at least, might not be mutually related in so simple way. But, 
nevertheless, the fact that the singularity spectra for Tg are wider than their 
counterparts for Gs and shifted towards larger a's is still in the spirit of the 
above-mentioned relation between the signed fx and unsigned fy processes. 

Our results suggest that as regards the original time series of price increments 
Ps{ts{i)) without any unfolding of the time axis, in that case we can well deal 
with a fully two-dimensional multifractal process, multiscaling both in the 
argument ts{i) and in the value Ps{ts{i))- This may suggest that both the price 
increments and the waiting times contribute to the multifractal properties of 
the time deformation 6{t) and of the fixed- At log- returns. 

Finally, Figure 9 shows a comparison of the spectra f{a) for the original Tg 
(solid), the reshuffled (dot-dashed) and the surrogate time series (dashed). 
Contrary to Fig. 5, here the randomized data displays completely different 
behaviour, being not only narrowed but also systematically shifted towards 
a = 0.5. This is not surprising, however: reshuffling removes the strong auto- 
correlation of the inter-trade time intervals completely. It can also be noted 
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Fig. 9. Comparison of the original and randomized series of waiting times: original 
(solid), reshuffled (dashed) and surrogate (dotted) time series. Multifractal spectra 
T{q) (left column) and singularity spectra f{a) (right column) are presented for 
three different stocks. Note the maximum positions for the surrogates with the 
linear correlations being preserved. 

that the spectra for the surrogate signals preserve the positions of the max- 
ima at 0.8; this is a natural consequence of the strong linear correlations still 
present in the surrogates. 



4 Conclusions 



We study the multifractal properties of the most liquid stocks from the Ger- 
man stock market. Our original data consisting of the recordings of time and 
price at which all the transactions took place, allowed us to separate the com- 
plete process of the stock trading into its pure price and pure time components. 
Since both components can contribute to the fixed-At volatility and returns, 
studying the properties of these components can help to identify the micro- 
scopic sources of the observed multifractality of the fixed-delay returns. We 
show that both the signals for the transaction-to-transaction price increments 
and for the inter-trade waiting times exhibit the characteristics that can be in- 
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terpreted in terms of multifractality. Its degree expressed by e.g. the widths of 
the singularity spectra /(a) varies across different stocks but these properties 
are entirely company-specific and are not related to industry sectors, com- 
pany size, average transaction frequency or any other characteristics of this 
type. The multifractal properties of Gg and Tg are of different nature; though 
on a qualitative level the corresponding spectra can be related. The relevant 
relation [16] between the /(a) spectra for the returns and for the multifrac- 
tal time 9{t) does not however apply fully quantitatively. This is because the 
c.d.f. of the microscopic price increments and the c.d.f. of the inter-trade time 
intervals are not related through the Brownian motion as required in [16]. If 
the price and the time components of the trading are independent or they at 
most weakly depend on each other, the compelling next step in this kind of 
analysis is to perform a fully 2-D multifractal approach in order to link the 
fractal nature of the price and the time increments into one unified frame. 
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